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ABSTRACT 

We investigate blending, binarity and photometric biases in crowded-field CCD imaging. 
For this, we consider random blend losses, which correspond to the total number of stars left 
undetected in unresolved blends. We present a simple formula to estimate blend losses, which 
can be converted to apparent magnitude biases using the luminosity function of the analyzed 
sample. Because of the used assumptions, our results give lower limits of the total bias and 
we show that in some cases even these limits point toward significant limitations in measuring 
apparent brightnesses of "standard candle" stars, thus distances to nearby galaxies. A special 
application is presented for the OGLE-II BVI maps of the Large Magellanic Cloud. We find a 
previously neglected systematic bias up to 0™2-0'.^3 for faint stars (V ^ 18™0 — 19™0) in the 
OGLE-II sample, which affects LMC distance measurements using RR Lyrae and red clump 
stars. We also consider the effects of intrinsic stellar correlations, i.e. binarity, via calculating 
two-point correlation functions for stellar fields around seven recently exploded classical no- 
vae. In two cases, for V1494 Aql and V705 Cas, the reported close optical companions seem 
to be physically correlated with the cataclysmic systems. Finally, we find significant blend fre- 
quencies up to 50-60% in the samples of wide-field exoplanetary surveys, which suggests that 
blending calculations are highly advisable to be included into the regular reduction procedure. 

Key words: techniques: photometric - methods: statistical - binaries: eclipsing - binaries: 
visual - stars: oscillations - stars: planetary systems - stars: novae, cataclysmic variables 



1 INTRODUCTION 

Classical crowded field photometry attempts to detect and mea- 
sure brightnesses of individual stars that are heavily affected by 
the presence of close neighbours. For ground-based observations, 
crowding depends on the angular density of objects and the at- 
mospheric seeing conditions. Stellar blending (unresolved imaging 
of overlapping stars) can be a significant component of the total 
ambiguity known as the confusion noise ~ see Takeuchi & Ishii 
(2004) for an excellent historic review and a general formulation 
of the source confusion statistics. Recent studies for which blend- 
ing was important include measuring: the luminosity function of 
individually undetectable faint stars (Snel 1998); the extragalactic 
Cepheid period-luminosity relation (Mochejska et al. 2000, Fer- 
rarese et al. 2000, Gibson et al. 2000); and Cepheid light curve 
parameters (Antonello 2002). Extensive investigations can also be 
found about blending and microlensing surveys (Alard 1997; Woz- 
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niak & Paczynski 1997; Han 1997, 1998; Alcock et al. 2001). Stel- 
lar blending in general is difficult to model, because significant con- 
tribution may be due to physical companions, which are common 
among young stars, including Cepheids (Harris & Zaritsky 1999, 
Mochejska et al. 2000). 

Here we attempt to determine the effects of random blending 
with a new approach that includes corrections for the excess of dou- 
ble stars. This work was motivated by recent cases in which stellar 
blending played a degrading role. For instance, wide-field photo- 
metric surveys of the galactic field are characterized by confusion 
radii of 10-20" (Brown 2003) and can suffer from strong blend- 
ing, even in regions far from the galactic plane. Other examples 
include the presence of a close optical companion of the classical 
nova V1494 Aql (at a separation of about l'.'5), which heavily af- 
fected late light curves of the eclipsing system (Kiss et al. 2004). 

Inspired by these problems and the availability of deep, all- 
sky star catalogues like the USNO Bl.O (Monet et al. 2003), we 
decided to perform simple calculations in terms of observational 
parameters such as the confusion radius - constrained by the see- 
ing or the pixel size of the detector - and stellar angular density. To 
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Table 1. Biases in apparent magnitude (rriobs) amplitude of variation 
(Aobs) for blending stars of magnitude difference Am. 
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the deep OGLE-II BVI map of the Large Magellanic Cloud and 
galactic novae in the USNO Bl.O catalogue. Case 2 corresponds 
to the exoplanetary surveys with very large confusion radii (up to 
20"), which can be well characterized using whole-sky star cata- 
logues of a magnitude better angular resolution. Case 2 is much 
simpler because one can always check whether an interesting object 
is a blend or not. However, if ensemble properties of large groups 
of stars are considered, blending must be taken into account be- 
cause high fractions of stars are in fact blends when observed at 
large confusion radius. 



emphasize the importance of the problem, we present magnitude 
and amplitude biases for unresolved blends of Am=0, 1, 2 and 3 
mag in Table Q Though fairly trivial, the numbers clearly show that 
even for Am=3 mag, blending can affect brightness and variabil- 
ity information to a highly significant extent. Also, these numbers 
span the magnitude ranges in which we are interested: usually 3 
to 5 mag wide samples are considered for the chance of random 
blending. This is a different approach compared to the case of gen- 
eral blending, which includes very faint blends, too. For instance, 
Han (1997) showed that for a typical I'.'S seeing disk towards the 
galactic Bulge, model luminosity functions of the Milky Way pre- 
dict ~36 stars within that area; of these, only 0.75% is expected 
to be brighter than 18 mag. Obviously, those faint blends do not af- 
fect photometry of bright stars. This study focuses on a much more 
specific problem: for a given field of view and confusion radius, 
what can be derived about random blending probabilities from the 
observed stellar angular distributions? Can we assign systematic 
errors based on these probabilities? 

The paper is organized as follows. In Sect. 2 we discuss blend- 
ing in random stellar fields. We present a simple formula to es- 
timate random blend losses, which has been tested by extensive 
simulations. We also discuss the effects of visual double stars by 
determining the two-point correlation. As an application, in Sect. 3 
we investigate probable photometric biases in deep OGLE-II BVI 
observations of the Large Magellanic Cloud (Udalski et al. 2000). 
In Sect. 4 we analyse stellar fields around seven recently erupted 
classical novae, all located in densely populated regions near to the 
galactic plane. We investigate blending rates in photometric survey 
programs HAT (Bakos et al. 2002, 2004), STARE^ (Brown 2003, 
Alonso et al. 2003) and AS AS (Pojmanski 2002) in Sect. 5. Con- 
cluding remarks are given in Sect. 6. 



2 RANDOM BLENDING AND BINARITY: BASIC 
RELATIONS 

Hereafter we distinguish two cases that need different approaches: 

Case 1: There is only one dataset to analyse for blending prob- 
ability, characterized by a typical confusion radius; no additional 
information exists based on a catalogue or high-resolution imaging 
with much smaller confusion radius. 

Case 2: We can compare the observations with an additional source 
of data, which can be considered as unbiased by random blending. 

Case I refers to those studies where the completeness of a cata- 
logue is investigated or where the observations were deeper than 
any existing catalogue. Examples discussed in this paper include 
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2.1 Random blend loss 

The idea behind our calculations is the following. When taking an 
image of a very crowded field, the detection efficiency is limited 
by the confusion radius (rc), which is the smallest angular distance 
between two resolvable stars. If the distance between two neigh- 
bours is smaller than Vc, we detect only one object. This means the 
number of detected stars, Nd, will be smaller than A'^ by the num- 
ber of objects in an area SS — ivr'^. The difference A'' — is 
the blend loss. Because of this, the detected "stars" will appear to 
be brighter than they really are. Our aim is to estimate the number 
of stars lost due to random merging and the corresponding mean 
magnitude biases. 

Let us consider a sample of A'^ single stars spread randomly 
over a field with area S. The mean number of stars within an area 
SS is 

5N = N ^ = nSS, (1) 

where n is the number density. If we assume that blend losses are 
predominantly due to close pairs of stars then the mean number of 
close neighbours within 5S can be expressed as ^nSS. In other 
words, we lose on average ^nSS stars for every SS area element 
that has been found to contain a star. Consequently, for a detected 
set of Nd stars, the total blend loss will be 

N-Nd = ^NdN^ (2) 

Rearranging this and introducing x = SS/S give the estimated 
total number of objects: 

Nd 



N 



l~^NdX 



(3) 



Note that - can be thought of as the number of resolution elements 
in the image. 

2.2 Monte Carlo simulations 

Equation 3 allows us to estimate the actual number of objects in 
our image, based on the number we have detected and the size of 
our resolution elements. We tested Eq. 3 with Monte Carlo simula- 
tions. Two different survey areas were chosen to mimic real obser- 
vations: 0.1 square degree with 1,000 objects and 1 square degree 
with 10,000 objects. This way the surface density was kept constant 
and the effects of statistics could be checked. The confusion radius 
was varied between I" and 40". For each simulation, we filled the 
survey area with randomly placed points and retained only those, 
that fell outside the confusion area for all previously placed points 
(i.e. after placing the first point, the second one was kept only if it 
was outside the confusion area of the first point; the third one was 
rc/ijeptetii^tiiilava&oill side the confusion areas of the first two points; 
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and so on). Every simulation was repeated one hiundred times and 
the average blending losses were compared to those of calculated 
via Eq. 3. 

The results are shown in Fig.Q Fractional losses range be- 
tween 0.1% and over 65%, with practically no difference between 
1,000 and 10,000 objects. The calculated blend losses are in ex- 
cellent agreement with the true values, except for the largest con- 
fusion radii. We stopped the simulations when the integrated con- 
fusion area reached the full survey area, because after that point, 
one can place infinitely many objects without detecting them. This 
is why the loss tends to be underestimated after r = 30". Also, 
this is where the mean number of objects within a confusion area 
exceeds 2, so that triplets and higher multiplets are no longer neg- 
ligible. Keeping these restrictions in mind, it is remarkable how 
well the blend losses agree with Eq. 3. However, one question has 
to be addressed before we turn to real datasets, for which binarity 
may be significant: to what extent can we assume that seemingly 
random stellar fields can indeed be approximated by a uniform ran- 
dom Poisson process? The answer can be found by checking the 
distance distributions of all pairs in a real sample. 

2.3 The effects of double stars 

We compare two histograms in Fig. |2| The thick black 
line shows the pair-distance distribution for USNO Bl.O stars 
around Nova (V705) Cas 1993 (21,520 stars with 14™00<red 
mag<17™99), while the grey line is for a simulation with 21,600 
random points with the same field of view. As expected, the over- 
all shapes of the distributions closely follow a Poisson-distribution, 
with slight boundary effect for the large distances (comparable to 
the diameter of the survey field), causing a somewhat sharper cut- 
off than a pure Poissonian. 

We have three conclusions based on these kinds of compar- 
isons for fields in the Milky Way: 

1 . The overall distributions agree very well with the pure random 
simulations. There are slight indications for different shapes (see, 
e.g., the upper panel of Fig.|2|around r = 1500" and 3000"), but 
the differences hardly exceed the intrinsic scatter of the data. 

2. The main disagreement occurs for the smallest pair- 
separations (typically for r < 10"), which must be due to a large 
number of binary stars. For example, in the field around V705 Cas, 
the pair excess suggests about 5000 visual double stars with sep- 
arations less than 7". In denser regions even higher numbers can 
occur: the 1 deg^ USNO Bl.O field around Nova (V4745) Sgr 2003 
contains about 20,000 double stars with red mags<17™0 and sep- 
arations under 10". Considering that the Washington Double Star 
Catalogue (Mason et al. 2001^) contains approximately 100,000 as- 
trometric double and multiple stars, it is evident that most faint bi- 
nary stars have yet to be identified and catalogued (see also Nichol- 
son 2002). 

3. The contamination by physical binary stars can therefore be a 
significant factor that may distort the results of random blending. 

Star-star correlations can be taken into account with the two- 
point correlation function, S,{r), which is defined by the joint prob- 
ability of finding an object in both of the surface elements SSi and 
5S2 at separation ri2: 

P = n^ 5 Si 5 S2ll + ar 12)]. (4) 

^ A regularly updated version can be found at 
|http : / /ad ■ usno ■ navy .mil /wds/| 
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Figure 1. Blend losses as a function of confusion radius in 0.1 deg^ (top) 
and 1 deg^ (bottom) fields, where the stellar density was set to 10'* star/deg^ 
in both fields. The heavily overlapping black dots represent one hundred 
simulations at each radius, while open circles connected with the solid lines 
show the predicted losses from Eq. 3. 



Note that homogeneity and isotropy are assumed, so that ^ is a 
function of the separation alone (Peebles 1980). If we can calcu- 
late this excess probability for pairs, then multiplying random blend 
losses from Eq. 3 by (1 -I- C('"c)) will correct the blending rate for 
the double stars. Neglecting higher correlation functions is a rea- 
sonable assumption, because the majority of stars belongs to single 
or double systems (about 90%, Abt 1983). Also, it is clear that even 
with this correction, the calculated blending rate will be a lower 
limit, because we do not have any information on the frequency of 
unresolved close binaries. 

We can determine ^(r) using the recipe that leads to Eq. 
(47.14) in Peebles (1980): place Nt points at random in the sur- 
vey area; let np{t) be the number of pairs among these trial points 
at separation r to r + 5r, and let Up be the corresponding number 
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Figure 2. Pair-separation distributions for a 1 deg^ field around V705 Cas 
(black line) and a random simulation (grey line). The lower panel is a zoom 
of the upper one showing the presence of about 5,000 resolved double stars. 



of pairs in the real catalogue of A'^ objects. Since ^ = for the trial 
points, the estimate of ^ for the data is: 



1 + 



(5) 



We show an example in Fig.|3| which in our case is simply 
the ratio of the two histograms plotted in Fig. |2| (because A'^t was 
equal to A^). The statistical uncertainty of the first point at r = 1" 
is ±0.01, far less than the symbol size, but a much larger system- 
atic error is likely because that radius is comparable to the image 
scale of plates on which the catalogue is based. The dashed line in 
Fig.|3]shows a power-law fit (^(r) ~ r") for the linear regime of 
the log-log plot. This diagram suggests that for the given dataset of 
21,500 stars, the probability of finding two stars with 1 arcsec sep- 
aration is about 20 times larger than for the pure random case. For 
a 1" confusion radius Eq. 3 gives A'^ — Nd = 56, which means the 
corrected blending loss, including the measurable fraction of dou- 
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Figure 3. The two-point congelation function for the 1 deg^ field around 
V705 Cas. The dashed line shows a power-law fit to the data between 2 and 
8 arcsecs (1 + ^(r) ~ r"-*-'^). 



ble stars, is about 1 120 stars. In other words, the chance of being a 
blend in this dataset is at least 5.2%. 

To conclude, it is always advisable to check star-star correla- 
tions before applying Eq. 3, and only if ^(r) is close to zero can we 
assume that double stars do not contribute significantly to blending 
(in practice, it means that the observations did not resolve binaries, 
probably because of the large distance). 



3 MAGNITUDE BIAS IN THE OGLE-II OBSERVATIONS 
OF THE LARGE MAGELLANIC CLOUD 

Our first application is an analysis of the OGLE-II BVI observa- 
tions of the Large Magellanic Cloud (Udalski et al. 2000). These 
data were obtained with the 1.3m Warsaw telescope at the Las 
Campanas Observatory for more than 7 million stars in the cen- 
tral 4.5 deg^ of the LMC. The completeness of the resulting cata- 
logue is high down to stars as faint as _B ~ 20 mag, V k, 2Q mag 
and I ~ 19.5 mag. The median seeing of the observations was 
I'.' 3, with no observations made when the seeing exceeded l'.'6- 
I'.'S (Udalski et al. 2000). Here we estimate blend losses in typical 
fields in the LMC and determine the corresponding magnitude bi- 
ases in V and /. 

This is a Case 1 blending, since no existing catalogue has a 
smaller confusion radius than the OGLE-II observations, except 
small field-of-view observations with the Hubble Space Telescope 
(Olsen 1999), which were used to estimate MACHO blend biases 
by Alcock et al. (2001). To investigate possible photometric bi- 
ases, we downloaded a representative field from the OGLE public 
archived We chose LMC_SC5, centered at RA(2000)=5''23™48', 
Dec(2000)=-69°41'05", which contains about 460,000 stars. The 
majority of stars are fainter than 18 mag (both in V and /), so 
we restricted the sample to stars between 18 mag and 21 mag. We 
used the two-point correlation function to constrain possible bina- 
rity effects and also to estimate the confusion radius. Because the 
stellar density varies gradually across the bar of the LMC, we se- 
lected several subsamples within which the density was constant. 
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Figure 4. Two-point correlation functions in V and / in Region 1 . 



Table 2. Star counts, estimated total numbers of stars, blend losses and 
blending frequencies for two OGLE fields in the LMC. 



Region 


band 


Nd 


N 


N-Na 


{N - Nd)/Na 


1 


V 


7598 


9151 


1553 


20% 


2 


V 


5354 


6081 


727 


14% 


1 


I 


7728 


9340 


1612 


21% 


2 


I 


6073 


7026 


953 


16% 



Here we discuss two of them, one in the middle of the bar (here- 
after Region 1, 80.7 < RA[°] < 80.8, -69.7 < Dec[°] < -69.6) 
and one further south (hereafter Region 2, 80.7 < RA[°] < 80.8 
and —70 < Dec[°] < —69.9). Other subsamples yield practically 
identical characteristics for blending systematics. 

In Fig.|4]we show the two-point correlation functions for Re- 
gion 1 in y and I bands. The step size in radius was 5r = 0'.'2, and 
we plotted the functions with linear scaling, so that changes around 
the confusion radius could be noticed more easily. We do not see 
a significant excess amount of correlated pairs, which is not sur- 
prising given the distance to the LMC. There is a slight rise of the 
correlation function for l'.'4 < r < 3" that may be attributed to the 
widest binary stars. However, the correlation goes down quickly for 
r < I'.'O. Since the best seeing of the OGLE images reached O'.'S 
(Udalski et al. 1998), we adopt this value as the confusion radius. 
Our choice is also consistent with "critical radius" chosen as 0'.'75 
by the OGLE-team, within which objects were treated as identical 
(Udalski et al. 1998). The blend losses estimated using Eq. 3 are 
shown in Table|2| 

Apparently, about one in five-six stars is a blend having an 
unresolved companion that is within 3 mags in brightness (since our 
data sample covers the range 18 to 21 mag). We have to stress that 
these numbers are only lower limits because of the incompleteness 
of the OGLE data below ~20 mag, which means the systematic 
errors we derive are still underestimated. Our main purpose is to 
point out the existence of surprisingly large systematic errors that 
can be determined from a single dataset alone. Deeper data (like 



in Alcock et al. 2001) would only shift these systematics toward 
larger values, but that is beyond the scope of the present paper. 

Having estimated the fraction of blended stars, we are now in a 
position to calculate the extent to which they will bias the photome- 
try. We assume that the probability of a star blending with a neigh- 
bour does not depend on its brightness, so that blended and un- 
blended stars have the same luminosity function (LF) within each 
region. This is likely to introduce a slight systematic error, but we 
estimate this to have a marginal effect on the final outcome of the 
investigation. Upper panels in Fig.|5|contain the normalized V and 
I luminosity functions (^^g LF(m) = 1), where the prominent 
peaks at V ~ 19 mag and / « 18 mag correspond to the red 
clump (Udalski 2000). The /-band decline for J > 20 mag shows 
the decreasing completeness (see Table 4 in Udalski et al. 2000), 
which is less severe in V. The next step was to carry out a Monte 
Carlo simulation for each region and each filter: we placed N stars 
at random into the given survey area, where A'' is the corresponding 
number in the fourth column of Table |2| The luminosity function 
was set to the one shown in Fig. |5| Then we determined which 
stars were "blends" in the random set and calculated the integrated 
magnitude of each blend. The magnitude difference between the 
blend and the brightest star within it was assigned as a bias at the 
integrated magnitude. Finally, these bias values were averaged for 
every 0. 1 mag bin of the luminosity function (unblended stars were 
taken into account by adding zero bias). The whole procedure was 
repeated one hundred times and the results are shown in the lower 
panels of Fig. |5| Small changes to the input luminosity functions 
affected the final bias distributions at a level lower than the scatter 
visible in the plots. 

We can draw several conclusions based on Fig.|5| Firstly and 
most importantly, blending introduces systematic errors as large as 
CK°2-0^3 (for instance, the mean bias between V = 18 — 19 mag 
is 0'"25 in Region 1 and 0^ 18 in Region 2), which can have serious 
consequences if neglected. Details of OGLE-II reduction, tests of 
photometric accuracy and incompleteness were published in Udal- 
ski et al. (1998) and they did not mention random blending as a 
possible source of systematic biases. Very recently, Alcock et al. 
(2004) presented an analysis of first-overtone RR Lyrae stars in 
the MACHO database, for which systematic magnitude biases due 
to blending were also estimated. At ~ 19^3, where RR Lyrae 
stars concentrate, Alcock et al. (2004) arrived at AV ranging from 
— CPllto— CK"21 for various assumptions (they made artificial star 
tests for various stellar densities and calculated the differences be- 
tween the input and recovered magnitudes of artificial stars). Our 
results are based on a very different approach applied to a very 
different dataset, and the agreement suggests that random blend- 
ing must be taken into account in such high-level crowding as 
the present one. Secondly, the gradual decrease of the bias to- 
wards fainter magnitudes shows the effect of incompleteness: as 
was shown by Alcock et al. (2001) for the MACHO database, sys- 
tematic errors increase monotonically over the examined magni- 
tude range. For that reason, the calculated biases between 18 and 
19 mag can be considered as lower limits for the fainter magni- 
tudes; the magnitude range of RR Lyrae stars is therefore very 
likely to suffer from a bias up to (f!^2-G^3, which has been ne- 
glected in the past. This leads to our third conclusion: results based 
on the OGLE-Il data that supported the LMC "short" distance scale 
((m — A/)lmc ~ 18™3) are likely to suffer from this systematic 
error. For example, Udalski (2000) calibrated the zero-point of the 
distance scale using RR Lyrae variables and red clump stars to be 
(m - M)lmc = 18™24 with 07^07 systematic error. However, 
blending correction significantly decreases the apparent conflict be- 
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Figure 5. The luminosity functions in Region 1 and 2 (upper panels) and tlie calculated systematic magnitude biases (lower panels). 



Table 3. Seven bright nova outbursts with reported progenitor candidates in 
the last decade. 



Star (year) 


Progenitor details 


Ref. 


V1494Aql (1999/2) 


Red mag 15™ 6 


1 


V2275 Cyg (2001) 


Red mag 18™ 8 (USNO A2.0) 


2 


V382 Vel (1999) 


Quiescent magnitude V = 16™56 


3 


V4743 Sgr (2002/2) 


Red mag 16™ 7 (USNO A2.0) 


4 


V4745 Sgr (2003) 


DSS red plate, R ^ 17™9, apparendy 






blended with a faint companion to SE 


5 


V705 Cas(1993) 


Candidate precursor of red mag about 
17™0, northern component of a merged 






pair with separation about 2 arcsec 


6, 7 


V723 Cas(1995) 


DSS red plate, mag 18-19 


8 



References: 1 - Pereira et al. (1999); 2 - Schmeer (2001); 3 - Platais et al. 
(2000); 4 - Haseda et al. (2002); 5 - Brown et al. (2003); 6 - Skiff et al. 
(1993); 7 - Munari et al. (1994); 8 - Hirosawa et al. (1995) 



tween this and the "standard" value (18™5, e.g. Alves 2004), which 
shows the importance of blending biases. 

It must be stressed that ground-based observations of more 
distant galaxies may suffer from even larger systematic errors. De- 
tecting variable stars by the image subtraction method circumvents 
the problem (see a recent application by Bonanos & Stanek 2003), 
but measuring apparent magnitudes (and hence distance) requires 
use of space telescopes. Calculations like the present one can help 
estimate systematics but more reliable results for nearby galaxies 
need high-resolution observations. 



4 FIELD STAR DISTRIBUTIONS AROUND SEVEN 
RECENT NOVAE 

Our next case study was inspired by the independent discovery of 
the close optical companion of nova V1494 Aql (Kiss et al. 2004). 
We selected several bright nova outbursts in the last decade for 
which progenitor star candidates have been identified in Digitized 
Sky Survey plates or deep star catalogues (most often in the USNO 
Al.O and A2.0 catalogue releases). The available details on these 



candidates are summarized in Table|3| We examined the following 
questions: How probable is it that the candidates are unrelated stars 
located only by chance at the nova coordinates? In cases where a 
star was found a few arcseconds from the nova position, (like for 
V1494 Aql and V705 Cas), how probable is is that they are physi- 
cally related to the nova system? 

To find out the answers, we downloaded and analysed 1 deg^ 
USNO Bl.O fields around each nova using the USNO Flagstaff 
Station Integrated Image and Catalogue Archive Service These 
fields provide a complete coverage down toV = 21™0, 0'.'2 astro- 
metric accuracy at J2000, 0.3 mag photometric accuracy in up to 
five colours and 85% accuracy for distinguishing stars from non- 
stellar objects (see Monet et al. 2003 for more details). The latter 
issue is less relevant in these fields because of the low galactic lat- 
itudes. In two cases we found incompleteness of the data, either 
in the form of large empty areas in the given field (V4745 Sgr) or 
sudden jumps in the stellar density within certain rectangular areas 
(V4743 Sgr). Both fields are located in the densest regions of the 
Milky Way, which explains the difficulties of measuring wide-field 
Schmidt-plates in such crowded fields. We took these into account 
in our analysis (e.g. in the latter case we omitted patches of lower 
density, yielding a decreased survey area). 

Since this is Case 1 blending, we calculated random blend 
losses for each nova by choosing all stars within 1-1.5 mag in 
brightness of the reported progenitor (when the accuracy was 
worse, we chose the wider range). For this, we used red magnitude 
data in the catalogue. We assumed a l'.'5 confusion radius, based 
on the pixel scale of the scanned Schmidt-plates (usually about 
O'.'9/pixel, Monet et al. 2003) and the fact that the optical com- 
panion of VI494 Aql located at l'.'4 was unresolved in the data. 
Two-point correlation functions (Fig. |6} were determined via Eq. 
5, filling up the survey areas by 25,000 points at random. Blend 
losses were then multiplied by the two-point correlation functions 
at r = l'.'5. Finally, we calculated the chance of finding a random 
blend within l'.'5 of the nova coordinates. We summarize the re- 
sults in Tableland in Fig.|6|(for V705 Cas we got essentially the 
same as in Fig.|3|for a wider magnitude range, so that there was no 
reason to repeat the plot). 

^ |http : / / www ■ nof s ■ navy .mil /data/ fchpix/| 
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Figure 6. Two-point correlation functions around six novae. 



Table 4. Statistical properties of the examined datasets. Abbreviations: A^^ - number of stars within the specified magnitude range; Nijn < m^) - total 
number of stars blighter than the fainter limit; AAf - random blend loss via Eq. 3; ^C^c) - the two-point coiTelation function at = 1''5; blend rate - 
blending occun'ence rate, corrected by 1 + ^(r); Al - the probability that the pre-outburst star coincides with the nova by chance; A2 - the probability that 
the optical companion is correlated with the nova itself. 



Star 


Mag range 




N{m < rriu) 




1 + 


blend 
rate 


Al 


A2 


VI 494 Aql 


15-16 


8,333 


13,229 


19 


4 


0.9% 


0.5% 


75% 


V2275 Cyg 


18-19 


22,249 


47,255 


137 


1.5 


0.9% 


1.2% 




V382 Vel 


16-16.6 


6,125 


19,291 


10 


15 


2.4% 


0.3% 




V4743 Sgr 


16-17 


10,775 


20,165 


32 


5 


1.5% 


0.6% 




V4745 Sgr 


17-18.5 


35,997 


62,042 


485 


3 


4% 


2.7% 




V705 Cas 


16.5-17.5 


8,066 


17,374 


18 


10 


2.2% 


0.5% 


90% 


V723 Cas 


18-19.5 


18,932 


31,302 


99 


3 


1.6% 


1% 





It is evident from Fig. |6| that the samples are dominated by 
binary stars in every field for separations under ~ 10". The highest 
correlation was found around V382 Vel, where the number of 1- 
arcsec doubles exceeds the number expected for a random field by 
a factor of 23! In other words, at least 22 of every 23 1-arcsec pairs 
are physically related (binary) stars. It is also prominent that there 
is a slope decrease in every correlation function for the smallest 
radii, which we interpret as a result of confusion losses. For that 
reason, we estimated ^(rc = 1"5) (6th column in Table |4j from 
extrapolated linear fits over the log-log plot of the 2-6 arcsecs data, 
rather than from interpolations between 1 and 2 arcsecs. 

The numbers in Table |4] clearly show that random blend loss 
stays well below 1%; even for the worst case, V4745 Sgr, it is only 
1.4% of Nd- After multiplying by the estimated two-point correla- 
tion function values, the blending rate still remains around 1-2%. 
The eighth column in Table|4|(Al) confirms that progenitors can be 
identified with high confidence, even in the central regions of the 
Milky Way. 

The most interesting numbers can be found in the last column 
of Table |4| They follow from the probability interpretation of the 
two-point correlation function. Since (1 + reflects the excess 



probability for finding a pair in a sample at distance r, 1 + ^(r) = 4 
(first row in Table|4} means that there are four times more pairs than 
in a random sample, so that for any given pair the chance of being 
correlated is 75%. Therefore, these data suggest that it is quite pos- 
sible that the optical companions of V 1494 Aql and V705 Cas form 
wide hierarchical triple systems with the cataclysmic components. 
Interestingly, for VI 494 Aql this is also supported by our knowl- 
edge of the optical companion. In Kiss et al. (2004) we assigned late 
F-early G as an approximate spectral type of the companion. The 
corresponding absolute magnitude for a main sequence star is about 
Mv =-l-4™0, which agrees within the error bars with the quiescent 
absolute magnitude of the nova (Kiss et al. 2004). Hence, the close 
apparent magnitudes of the pre-nova and the companion suggest 
similar distances, consequently the possibility of physical correla- 
tion. For V705 Cas, we do not have similar supporting evidence, 
but it is an intriguing possibility that these optical components are 
in triple systems with the novae. If confirmed, their presence could 
be used, for example, to derive independent distances via spectro- 
scopic parallaxes, which in turn would allow one to calculate accu- 
rate absolute magnitudes of the nova eruptions. 

This case study shows that seeing-limited images are not sig- 
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Figure 7. Blending rates for similar projects at four different galactic latitudes (HAT: +30° ; STARE: +2° ; ASAS(l): -60° ; ASAS(2): -20° ). The confusion 
radius is about 20" for each project. 



nificantly affected by random blending, even in the galactic plane. 
We now turn to the case of images with much lower spatial resolu- 
tion. 



5 WIDE-FIELD PHOTOMETRIC SURVEYS 

In the recent years there has been an explosion in the number of 
small and medium-sized robotic telescopes that monitor selected 
regions of the sky with various instruments (e.g., Hessman (2004) 
listed 80 different projects). We selected three representative small 
instruments that are currently running, to investigate Case 2 blend- 
ing. They are: 

(i) The Hungarian Automated Telescope (HAT) which, in its 
present status, observes stars between 7 = 9 — 12 mag in selected 
fields of the northern sky. The data were reduced with PSF-fitting 
photometry (Bakos et al. 2004). We downloaded the 1 deg^ USNO 
Bl.O field in Hercules (RA(2000)=17''36'", Dec(2000)=37°30', 
6« 30°). 

(ii) The STellar Astrophysics & Research on Exoplanets 
(STARE) project, observing between _R = 9 — 12"!" 5 in var- 
ious fields (Alonso et al. 2003, Brown 2003). We downloaded 
the 1 deg^ USNO Bl.O field, centered at RA(2000)=20''06™, 
Dec(2000)=36°00' (b ^ 2°). 

(iii) The All Sky Automated Survey (ASAS), which covers 
over 1 million stars in the southern hemisphere between V = 
8™5 — 15 mag. The data are reduced with simple aperture pho- 
tometry (Pojmanski 2002, 2004). We selected two 1-deg^ fields 



at random at two different galactic latitudes: RA(2000)=00 19™, 
Dec(2000)=-56°00' (b f» -60°) and RA(2000)=06''10'", 
Dec(2000)=-23°00' {b ^ -20°). 

Each project is characterized by rc ~ 20", so that their com- 
parison reveals the differences that depend on the galactic latitude. 
To match the photometric bands used in these projects, we exam- 
ined USNO Bl.O infrared magnitudes for the HAT-field and red 
magnitudes for the remaining regions. The main aim was to get 
an idea of the fraction of biased stars, because although crowding 
problems are well-known for these instrumental setups, we did not 
find any quantitative description of the issue. Therefore, we cal- 
culated blending rates for various m and Am values that covered 
the magnitude ranges of the projects. This rate was defined as the 
following ratio: 



blending rate 



p{m, Am, Tc) 



(6) 



where u{m, Am, rc) is the number of stars which have fainter 
neighbours within a distance and magnitude difference Am; 
i'(m) is the total number of m magnitude stars (in our case, de- 
fined as stars with apparent brightnesses within m and m + 1 
mag). We assumed negligible blend losses in the USNO Bl.O cat- 
alogue over the studied magnitude and separation ranges. The re- 
sults are shown in Fig.|7| We also determined two-point correlation 
functions, which showed that for these bright magnitudes the over- 
whelming majority of pairs closer than the confusion radius are 
physically related double stars. 
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Figure Q implies an alarming rate of blending, even for high 
galactic latitudes. We see that 10 to 20 percent of objects observed 
by the HAT and ASAS projects have blends within 3 mags, while 
up to 50% are affected near the galactic plane (STARE). Correlated 
pairs can make the situation quite bad even for the brightest stars, 
which means blending must be taken into account in every case. 
Presently available star catalogues offer a good opportunity to do 
that, thus we recommend to add blending information in all cases 
when finally reduced data are made accessible to the wider com- 
munity. Also, the implementation of the image subtraction method 
(Alard & Lupton 1998) is highly desirable in this type of project, 
because the photocenter of the variable source can be used to iden- 
tify which object within a blend is varying (Alard 1996; see also 
Hartman et al. 2004). 



6 CONCLUSIONS 

In this study we investigated the effects of random blending that in- 
volves stars of similar brightnesses, leading to significant biases in 
the measured apparent magnitudes and amplitudes. A simple for- 
mula (Eq. 3) was derived to calculate blend losses in a real cata- 
logue of stars based on the confusion radius, total number of de- 
tected stars and survey area. We showed that the two-point cor- 
relations must be included in calculations when studying galactic 
fields, where binary stars dominate for separations under 10". Out- 
side the Milky Way, we quickly lose the information on wide bina- 
ries, so that the pure random case applies. That is why the calcu- 
lated blending rate always puts only a lower limit to the full blend- 
ing. 

We discussed three different applications, which demonstrate 
the importance of these phenomena. The most interesting results 
were presented for the OGLE-II data of the Large Magellanic 
Cloud. With extremely high stellar angular densities, these observa- 
tions are much more biased by random blending than those fields 
in the galactic plane. Even though no observations were carried 
out for seeing larger than l'.'6-l'.'8, we estimated that 15-20% of 
the sample are affected by blends that are within three magnitudes 
in apparent brightness. The resulting magnitude biases can reach 
(K"2-0?^3, depending on the angular position within the LMC. This 
allows us to reconcile the OGLE-based distance moduli of the LMC 
(fj, Ri 18™ 3) and the generally adopted "standard" one (/i ~ 18™ 5). 
Random blend calculations are therefore highly desirable in ev- 
ery case when significant blending is expected, especially in ex- 
tragalactic observations. 

The analysis of star fields around seven novae showed that 
progenitor identification is very secure, even in the extremely dense 
fields of the Bulge. Furthermore, statistical evidence points toward 
the existence of wide hierarchical triple systems, in which the third 
component lies at 1-2" (i.e. several thousand AU) from the no- 
vae. Thus, it is highly desirable to measure proper motions for 
V1494 Aql and V705 Cas and their optical companions to confirm 
or reject the physical correlations. 

In our last example we demonstrated high-blending rates in 
typical wide-field photometric surveys for exoplanet transits. In 
certain cases, up to 50-60% of stars can be affected by blending 
objects within 3 magnitudes. We recommend that cross-correlation 
with appropriate star catalogues should always be done as part of 
the regular data reduction, to flag every possibly problematic star. 

Our presented method for calculating blend losses is very sim- 
ple, thus may not be applicable in certain cases. For instance, when 
the stellar angular density shows strong gradient, like in a globu- 



lar cluster or outer parts of a resolved galaxy, the basic assump- 
tion of homogeneity and isotropy will be invalid. A possible so- 
lution in such cases is to introduce a position-dependent density 
n — p{x, y), where x and y are the image coordinates, and to inte- 
grate all equations over the inhomogeneous and anisotropic sample. 
In principle, the numeric implementation of this generalisation is 
not too difficult. Another way to treat these cases is to split the data 
into small segments within which n ~ const, and, after calculating 
blending biases, take averages over the segments. 

To summarize, the experience with real datasets suggests that 
it is always highly advisable to estimate random blend losses when 
the star density and/or the confusion radius have relatively large 
values - Eq. 3 or deep star catalogues can be used to characterize 
any specific example. Where applicable, the upgrade from aperture 
and PSF-photometry to the image subtraction method is desirable 
to reduce some of the blending effects. 
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